CURVATURE BASED METHOD FOR SELECTING FEATURES FROM AN 
ELCTROPHYSIOLOGIC SIGNALS FOR PURPOSE OF COMPLEX 
IDENTIFICATION AND CLASSIFICATION 
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Cross-Reference to Related Application(s) 

This application is a continuation of U.S. Patent Application Serial No. 
09/703,269, filed on October 31, 2000, the specification of which is incorporated 
herein by reference. 

Field of the Invention 

The present invention relates generally to cardiac rhythm management 
systems, and more particularly, it pertains to a system and method of classification and 
detection of cardiac signals. 



Background of the Invention 

Analysis of cardiac signals, which is routinely performed in 
electrocardiography is generally based on visual inspection to quantify or qualify wave 
morphology for the purpose of identifying and classifying abnormal pattems. Certain 

20 morphological characteristics of commonly recorded signals have high diagnostic 

value. The shape and inter arrival times of R-waves recorded in the electrocardiogram 
generally provide a wealth of information about the state of the heart. Accordingly, 
automated approaches for identifying and classifying abnormalities in signals such as 
cardiac signals have sought to use a signal's significant morphologic characteristics. 

25 However, given the wide diversity of possible shapes for cardiac signals, it is 

usually not possible for an automatic approach to identify significant characteristics 
that can be used for unambiguous classification. Rather, automated classification 
approaches generally compare the entire morphological shape of a signal with the 
shape of similar signals with known abnormalities but without particular regard to the 

30 specific characteristics that the signals contain. Alternatively, automated classification 
approaches restrict the automated examination only to those signals which are 
essentially normal and use detailed metrics (for example QRS width, QT interval or 
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ST segment amplitude) of the essentially noraial morphology for classifying 
abnormalities. 

Despite its importance in the analysis of biologic signals, the automated and 
accurate identification and quantification of the significant morphological 
5 characteristics (for example turns, peaks, knees, inflection points, and the like) in any 
cardiac signal (both abnormal as well as normal) is still in a developing stage. Existing 
methods have used the concept of sharpness (for example to detect R-waves) but have 
had limited success. This is due in part to the overly simplistic mathematical treatment 
this concept has received, as reflected in the rudimentary algorithms used for these 

10 measurements. Most of the current detection methods rely on three point 

interpolations to measure sharpness. The simplest and most commonly used methods 
for measuring peaks of R-waves are based upon Taylor-series approximations to 
estimate the second derivative of the sensed signal. This formula utilizes highly local 
information (the point at the peak and its two close neighbors) ignoring nearby points 

15 which may contribute to signal peak. Other popular approaches utilize less local data, 
such as the peak and two adjacent extrema. All of these methods, which rely on three- 
point estimates of sharpness, may produce inaccurate estimates, if waveforms are 
complex or are contaminated with noise. Thus, a need exists for automated 
identification and classification of peaks, knees, inflection points, and the like in 

20 sensed cardiac signals that takes into account wave scale and complexity that can yield 
a more accurate estimate of peaks for identifying and classifying abnormalities in 
cardiac signals. 

Summary of the Invention 

25 The present subject matter provides a curvature based method of selecting 

features firom electrophysiologic signals for purpose of complex identification and 
classification. According to one aspect of the present subject matter, this is 
accomplished by sensing a cardiac signal (sensing the cardiac signal includes sensing 
complexes continuously on a real-time basis) and computing curvatures on a sample 

30 point-by-sample point basis (Xi, X2, X3, Xi ) on the sensed cardiac signal on a 
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continuous basis. In one embodiment, the curvature at the sample points Xi, X2, X3, 

Xi are computed by fitting a cubic least square error curve (using N number of 

sample points) to the sensed cardiac signal In this embodiment, N is an odd number, 
and the sample point (where the curvature is computed) is at a mid point of the N 
5 number of sample points. 

Also on a continuous basis, features of significant interest are extracted from 
the computed curvatures. In some embodiments, this is accomplished by comparing 

the computed curvatures at the sample points Xi, X2, X3, Xi to a set of 

predetermined threshold values. In some embodiments, the features are extracted 

10 based on computing features such as a time when the feature occurs, an amplitude of 
the feature, and other similar features at each of the sample points and comparing them 
to a set of threshold values. 

Also on a continuous basis a set of features associated with a first complex in 
the sensed cardiac signal are identified and separated from the continuously computed 

15 and extracted features upon detecting a second subsequent complex. The second 
subsequent complex is a complex that is adjacent to the first complex and occurs 
substantially inunediately after the first complex. This process of identifying and 
separating extracted features repeats itself from one sensed complex to another 
subsequent sensed complex on a real time basis. One reason for identifying and 

20 separating the set of features associated with the first complex is to prevent the 

features associated with the first complex from mixing with the features associated the 
second subsequent complex. Separating the features associated with first complex aids 
in classifying the sensed first complex. 

Next, the process includes identifying a fiducial feature from the separated set 

25 of extracted features associated with the first complex and aligning the separated set of 
features with respect to the identified fiducial feature. In one embodiment, fiducial 
feature is identified based on comparing the times when each of the separated set of 
features occur with a time when a complex associated with the separated set of 
features is detected on the sensed cardiac signal, and selecting a feature from the 

30 separated set of features that is closest in time to the time when the complex associated 
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with the separated set of features was detected. One reason for using a time when a 

complex is detected (such as R wave) in identifying the fiducial feature, is because the 
detection of a complex in a sensed cardiac signal is generally more reliable and 
consistent. 

5 Next, the process includes aligning the separated set of features around the 

identified fiducial feature. Aligning the separated set of features around the identified 
fiducial feature aids in normalizing each of the separated set of feature around a datum 
such as the associated detected complex, and further aids in comparing the separated 
set of features with a set of predetermined templates. 

10 Next the process includes comparing the aligned set of features to a set of 

predetermined templates to classify the associated complex. In some embodiments, the 
predetermined templates are a set of identified complexes associated with known 
cardiac arrhythmias that would assist in comparing and classifying the extracted set of 
features. In some embodiments, a therapy is provided to the heart based on the 

15 outcome of the classification. The above described process repeats itself on a 

continuous basis for a real-time classification of complexes fi*om the sensed cardiac 
signal. 

Other aspects of the invention will be apparent on reading the following 
detailed description of the invention and viewing the drawings that form a part thereof 

20 



Brief Description of the Drawings 
Figure 1 is a cardiac signal according to one aspect of the present invention. 
Figure 2 is a curvature figure derived fi-om the cardiac signal of Figure 1 
25 according to one aspect of the present invention. 

Figures 3A and 3B illustrate separating and aligning features derived fi^om the 
curvature figure of Figure 2 according to one aspect of the present invention. 

Figure 4 is a flow diagram illustrating generally one embodiment of operation 
of the present subject matter. 
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Figure 5 is a schematic/block diagram illustrating generally one embodiment 
of portions of a cardiac rhythm management system of the present subject matter. 

5 Detailed Description 

In the following detailed description, reference is made to the accompanying 
drawings which form a part hereof, and in which is shown by way of illustration 
specific embodiments in which the invention may be practiced. These embodiments 
are described in sufficient detail to enable those skilled in the art to practice the 

10 invention, and it is to be understood that the embodiments may be combined, or that 
other embodiments may be utilized and that structural, logical and electrical changes 
may be made without departing from the spirit and scope of the present invention. 
The following detailed description is, therefore, not to be taken in a limiting sense, and 
the scope of the present invention is defined by the appended claims and their 

15 equivalents. 

Analysis of cardiac signals, which is routinely performed in 
electrocardiography is generally based on visual inspection to quantify or qualify wave 
morphology for the purpose of identifying and classifying abnormal patterns. Certain 
morphological characteristics of commonly recorded signals have high diagnostic 

20 value. The shape and inter arrival times of R-waves recorded in the electrocardiogram 
generally provide a wealth of information about the state of the heart. Accordingly, 
automated approaches for identifying and classifying abnormalities in signals such as 
cardiac signals have sought to use a signal's significant morphologic characteristics. 
However, given the wide diversity of possible shapes for cardiac signals, it is 

25 usually not possible for an automatic approach to identify significant characteristics 
that can be used for unambiguous classification. Rather, automated classification 
approaches generally compare the entire morphological shape of a signal with the 
shape of similar signals with known abnormalities but without particular regard to the 
specific characteristics that the signals contain. Alternatively, automated classification 

30 approaches restrict the automated examination only to those signals which are 
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essentially normal and use detailed metrics (for example QRS width, QT interval or 
ST segment amplitude) of the essentially normal morphology for classifying 
abnormalities. 

Despite its importance in the analysis of biologic signals, the automated and 
5 accurate identification and quantification of the significant morphological 

characteristics (for example tums, peaks, knees, inflection points, and the like) in any 
cardiac signal (both abnormal as well as normal) is still in a developing stage. Existing 
methods have used the concept of sharpness (for example to detect R- waves) but have 
had limited success. This is due in part to the overly simplistic mathematical treatment 
10 this concept has received, as reflected in the rudimentary algorithms used for these 
measurements. Most of the current detection methods rely on three point 
interpolations to measure sharpness, may produce inaccurate estimates, if waveforms 
are complex or are contaminated with noise. 

15 General System Overview 

The present subject matter provides, among other things, a system of 
classifying a cardiac signal by computing a cxuvature versus time signal firom a sensed 
cardiac signal and comparing the computed curvature versus time signal with a set of 
predetermined templates. The present system also provides an automated 

20 identification and classification of peaks, knees, inflection points, and the like in 
cardiac signals that takes into account wave scale and complexity that can yield a 
more accurate classification of cardiac signals. 

The first step in the process includes computing continuously a curvature 
versus time signal from a sensed cardiac signal on a real-time basis. The process also 

25 includes extracting features (features are generally significant points of interest on the 
sensed cardiac signal) continuously on a real-time basis fi-om the computed curvature 
versus time signal. Next, the process includes continuously separating a set of features 
associated with a first complex upon detecting a second subsequent complex so that 
each individual complex can be identified and classified on a real-time basis. Then the 

30 process includes aligning and classifying the extracted and separated set of features by 
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comparing the separated set of features with a set of predetermined templates. Other 
aspects of the invention will be apparent on reading the following detailed description 
of the invention and viewing the drawings that form a part thereof 

Referring now to Figure 1, there is one embodiment of the first step in the 
5 process of sensing a cardiac signal 100 according to the present invention. Shown in 
Figure 1 are detection of heart beats associated with first and second complexes 110 
and 120 at times '^y' and 7^' associated with the sensed cardiac signal 100 on a real- 
time basis on a time line 130. Also, shown on the time line 130 are the various sample 

points Xi, X2, X3, Xi used in computing a curvature versus time signal on the 

10 first complex 1 10. In the example embodiment shown in Figure 1, the sample points 
can be some fixed interval of time, such as 2, 4, 6... milliseconds, at which curvatures 
are computed on the sensed cardiac signal. Also shown in this embodiment, is the 
detection of the complexes (1 10, 120, and so on) on a continuous basis on the sensed 
cardiac signal 100. 

15 Referring now to Figure 2, there is shown one embodiment of generating a 

curvature versus time signal 200 on a real-time basis derived fi-om the sensed cardiac 
signal 100 shown in Figure 1 according to the present invention. Figure 2 is also on 
the same time line 130 as in Figure 1. Curvature figure 200 shown in Figure 2 is a 
computed curvature (expressed in time) versus time signal generated continuously by 

20 computing curvatures at the sample points Xi, X2, X3, Xi on the sensed cardiac 

signal 100 shown in Figure 1. Curvature figures 210 and 220 are computed curvature 
versus time signals associated with first and second complexes 110 and 120 shown in 
Figure 1, respectively. In this example embodiment, the curvature versus time signal is 
computed by overcoming the dimensionality of the curvature by changing the 

25 dimensionality of the curvature to time. The reason for overcoming the dimensionality 
of the curvature is to convert the time versus voltage signal, which is time versus 
curvature, to a time versus time signal to facilitate easier computation of curvature (for 
example, when a curve exists in length-length space, the quantity has dimension of 
length and is traditionally called the radius of curvature, and if the curve exists in 

30 time-time space, the analogous quantity will be called the radius of time). In this time- 
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time space the curvature has a dimension of 1/Time. Also, shown in Figure 2 are some 
example sample points 230 at which curvatures are computed. 

Referring now to Figures 3 A and 3B, there is shown one example embodiment 
of extracting and aligning 300 sets of features associated with respective detected 
5 complexes. Figure 3A, illustrates continuously extracting features 320 from the 
computed curvature versus time signal 200 shown in Figure 2. Also shown in this 
example embodiment, are some sample times in milliseconds 330 when the extracted 
features occur on the time line 130. Figure 3 A is also on the same time as that of 
Figures 1 and 2. In the example embodiment shown in Figure 3 A, the sample times 

10 (where the curvatures are computed) have a fixed interval of 2 milliseconds. Also, in 
this example embodiment, the beat associated with first complex 1 10 is detected at 1 1 
milliseconds {ti) and the beat associated with the second subsequent complex 120 is 
detected at 101 milliseconds (/i^. 

Figure 3B shows separating sets of extracted features 310 associated with the 

15 complexes 1 10, 120 ... on a continuous basis from the continuously extracted features 
320 shown in Figure 3A. In the example embodiment shown in Figure 3B, a set of 
features 340 associated with first complex is separated upon detecting the second 
subsequent complex 120. Also shown are separating a set of features 350 associated 
with second complex upon detecting a third subsequent complex. This process of 

20 separating a set of features associated with a complex is done on a continuous basis. 
Features associated with first complex is separated from the features of subsequent 
second complex in order to (so that the features are not mingled with the extracted 
features of the second subsequent complex) classify the sensed first complex 1 10. The 
next step is to identify fiducial features 360, 370, and so on. In the example 

25 embodiment shown in Figure 3B, fiducial feature 360 is identified based on comparing 
times when the separated sets of extracted features associated with the first complex 
110 occur with a time when the heart beat associated with the first complex 1 10 is 
detected ('// ' ), and choosing a feature from the set of features 340 that is closest in 
time to the time when the heart beat associated with the first complex 1 10 is detected. 

30 In the example shown in Figure 1, the first complex 1 10 is detected at 1 1 miUiseconds. 
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The times when the separated set of features 340 associated with the first complex 110 
occur are 2, 4, 8, 10, 12, 14, 20, and 24 milliseconds as shown in Figure 3 A. Based on 
comparing the times when the separated set of features 340 occur with the time when 
first complex 1 10 was detected, the feature 360 that occurs at a time of 10 
5 milliseconds is the closest in time to the detected first complex 110 which occurs at 1 1 
milliseconds. Therefore, extracted feature 360 is selected as the fiducial feature fi-om 
the separated set of features 340 associated with the first complex 110. Similarly, the 
extracted feature 370 is selected as the fiducial feature for the separated set of features 
350 associated with subsequent second complex 120. This process repeats itself 

10 continuously on a real time basis. In some embodiments, the fiducial feature is 

identified fi-om the set of separated features associated with the first complex based on 
a predetermined deviation value. The predetermined deviation value is based on a 
sample point having an amplitude farthest firom a predetermined reference point. 

The next step is to align the separated set of features 340 associated with the 

15 first complex 110, and to classify the separated set of features 340 by comparing the 
aligned set of features with a set of templates. Separated sets of features 340 and 350 
are aligned to normalize the separated sets of feature 340 and 350 so that the aligned 
features aid in comparing the separated sets of features with a template to classify the 
sensed complexes 110 and 120. Figure 3B, illustrates one example embodiment of 

20 aligning the set of extracted features 340 (associated with first complex 110) with 
respect to an identified fiducial feature 360. In this example embodiment, the time 
when the identified fiducial feature occurs (at 10 milliseconds) is subtracted from each 
of the times when the separated set of features 340 associated with the first complex 
110 occurs (at 2, 4, 8, 10, ... milliseconds) to arrive at ahgned times -8, -6, -2, 0, 2, 4, 

25 10, and 14 milliseconds for the separated set of features 340. Similarly, the aUgned 
times -8, -2, 0, 2, 4, and 10 are computed for the separated set of features 350 
associated with the second complex 120. This process of aUgning separated set of 
features continues on a real time basis for the detected complexes on the sensed 
cardiac signal 100. It can be seen from the above illustrated computation, that the 

30 process aligning helps normalize the times associated with each of the separated sets 
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of features, and how it aids in comparing the separated set of features 340 and 350 
with a set of templates and in classifying the sensed complexes 110 and 120. 

Figure 4, illustrates a method 400 of identifying and classifying a sensed 
cardiac signal on a real-time basis. Method 400, as shown in Figure 4, begins with step 
5 410 by sensing a cardiac signal 410 from the one or more electrodes disposed in or 
around a heart 515 on a continuous basis. 

The next step 420 in the process comprises computing curvatures at sample 

points Xi, X2, X3, Xi on the sensed cardiac signal In some embodiments, the 

sample points Xi, X2, X3, Xi are fixed intervals of time at which the sensed 

10 cardiac signal is computed for its curvature. The sensed cardiac signal can be an atrial 
or ventricular signal of a heart. In one embodiment, the curvatures at sample point Xi, 

X2, X3, Xi is computed by fitting a cubic least square error curve to the sensed 

cardiac signal by using N sample points. In this embodiment, N is an odd number 
greater than or equal to 5. Also in this embodiment, the cardiac rhythm management 

15 system computes the curvature at a mid point of the N sample points. In one 

embodiment, the cardiac rhythm management system computes curvature (K) at a 
sample point Xi on the cardiac signal by using 5 sample points over which the cubic fit 
is (N being 5 in this example) made using 

20 K-(2C,/(1+Biy^) 
where 

Bi = T,Y(I-2) 4- T2Y(I-1) + T3Y(I) +.T4Y(I+1) + T5Y(R2) 
Ci = SiY(I-2) + S2Y(I-2) + S3Y(I) + S4Y(I-Hl) -h S5Y(I-h2) 
where Si, S2, S3 ... and Ti, T2, T3 ...are constant terms in the summation. In this 
25 embodiment, I corresponds to the sample point at which the curvature is being 

computed. For example, if the sample point is X3, then I will be equal to 3 and so on. 

In another embodiment, the system computes average curvature between two 
adjacent sample points by using a linear interpolation of the Bi and Ci between the 
two adjacent sample points (this is done by doing the integral of the curvature between 
30 the two adjacent sample points). This integration helps capture any curvatures that 
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become larger between the two adjacent sample points, so that the point-wise 
curvature values may more robustly represent the curvatures in the sensed cardiac 
signal. The system also integrates curvature versus time curve to find the area under 
the curve. This area under the curve represent the total change in angle between two 
5 points on the original cardiac signal. Thus, peaks, turns, inflections, and the like in the 
cardiac signals can be identified when the curvature versus time curve has peaks or 
nadirs. These peaks and nadirs denote points along the sensed cardiac signal where the 
tums are locally sharpest. Also, slow but steady turns on the sensed cardiac signal can 
be characterized by a small peak having a substantial area under the curvature curve. 
10 Thus, integrating curvatures between two adjacent sample points can help capture any 
curvatures that may become larger between any two adjacent sample points. 
In one embodiment, the cardiac rhythm management system computes 

curvatures at sample points Xi, X2, X3 Xj on the sensed cardiac signal on a 

continuous basis. 

15 Then the next step 430 in the process includes extracting features from the 

computed curvatures by comparing the computed curvatures with a set of 
predetermined threshold values. Figure 3A, illustrates extracting features fi:om the 

computed curvatures at the sample points Xi, X2, X3, Xi. In some 

embodiments, the set of predetermined threshold values are based on a previous 

20 curvature value, a first curvature value, and a curvature threshold limit. In one 

embodiment, the step of extracting features fi-om the curvature versus time signal is 
done on a continuous basis on the sensed cardiac signal as shown in Figure 3A. In 
some embodiments, the features are defined by one or more metrics. In this 
embodiment, the one or more metrics are defined as area under a computed curvature, 

25 a time of centroid of the area, and a value of original signal amplitude at a time of 
centroid of the area associated with the feature. 

Then the next step 440 comprises separating a set of features associated with a 
first complex upon detecting a second subsequent complex fi'om the sensed cardiac 
signal Separating of extracted features associated with the detected complex is 

30 illustrated in detail Figures 1, 2, 3 A and 3B. In some embodiments, the second 
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subsequent complex is detected based on continuously processing the sensed cardiac 
signal to produce a signal that is the absolute value of the first derivative of a sensed 
cardiac signal versus time. Then comparing the produced absolute value versus time 
signal with a predetermined decaying threshold value to detect an occurrence of the 
5 second subsequent complex. In some other embodiments, the second subsequent 
complex is detected when the absolute value of the first derivative exceeds a 
predetermined decaying threshold value. In some embodiments, the features 
associated with the first complex are stored in a memory until the second subsequent 
complex is detected, and at that point, the features associated with the first complex 

10 are replaced by the features associated with the second subsequent complex in the 
memory. The step of separating can comprise separating the set of features based on 
selecting features associated with the first complex having a predetermined time 
earlier than the detected second subsequent complex. In some embodiments, each of 
the set of features are further compared to an area threshold value to eliminate features 

15 less than the area threshold value. This can result in a set of revised features that can 
represent the first complex. 

Then the next step 450 in the process is to identify a fiducial feature firom the 
set of separated features associated with the first complex based on a predetermined 
deviation value. In some embodiments, fiducial feature is identified based on 

20 comparing the time when the first complex occurs with the times when the separated 
set of features associated with the first complex occur. The process of identifying a 
fiducial feature is discussed in more detail with respect to Figures 3A and 3B. 

Then the next step 460 in the process comprises aligning the set of separated 
features with respect to the identified fiducial feature. The process of aligning is also 

25 discussed in more detail in Figure 3B. The process of aligning helps normalize the 
times associated with each of the separated set of features and further facilitates in 
comparing the separated set of features with a set of predetermined templates to 
classify the detected complex. 

Then the next step 470 comprises comparing the aligned set of separated 

30 features with a set of predetermined templates. In one embodiment, the set of 
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predetermined templates are one or more known types of heart beat signals. In another 
embodiment, the set of predetermined templates consists of sets of template boxes or 
zones for each of the known types of heart beat signals. Generally, each template box 
or zone can have center amplitude, a time width, and an amplitude width. Also 
5 generally, the template boxes or zones for each known type of heart beat can consist of 
one or more template boxes or zones each of which can have different centers and 
widths. In this embodiment, the comparing step comprises assigning a score based on 
how well each feature matches with each of the set of template boxes or zones. In this 
embodiment, if the separated set of features associated with the first complex fails to 

10 match with the set of templates then the complex is assigned with an unknown heart 
beat classification. In one embodiment, the set of template boxes or zones are 
modified over time to reflect any changes in sets of separated features used in 
matching with the template boxes or zones. This can permit known beat type 
templates to follow slow trends or long-term changes in the cardiac signal. 

15 The next step 480 in the process comprises classifying the detected first 

complex based on the outcome of the comparison. Then the next step 490 comprises 
providing a therapy based on the outcome of the classification. In one embodiment, 
the next step includes guiding a therapy to a heart based on the outcome of the 
classification. In another embodiment, the next step in the process includes storing 

20 classifications for future diagnostic purposes. 

Figure 5, is a schematic/block diagram 500 illustrating one embodiment of 
portions of a cardiac rhythm management system 505, which is coupled to a heart 515. 
Cardiac rhythm management system 505 includes a power soxu:ce 580, a controller 
540, a sensing circuit 520, a therapy circuit 590, and a via node/bus 530. The sensing 

25 circuit 520 is coupled by a lead 5 1 0 to the heart 5 1 5 for receiving, sensing, and or 

detecting electrical heart signals. The sensing circuit 520 provides one or more sensed 
cardiac signals to the controller 540, via the node/bus 530. The controller 540 also 
controls the delivery of a therapy provided by the therapy circuit 590 and/or other 
circuits, as discussed below. 
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The controller 540 includes various modules, which are implemented either in 
hardware or as one or more sequences of steps carried out on a microprocessor or 
other controller. Such modules are illustrated separately for conceptual clarity; it is 
understood that the various modules of the controller 540 need not be separately 
5 embodied, but may be combined and/or otherwise implemented, such as in 
software/firmware. 

In general terms, the sensing circuit 520 senses a cardiac signal on a 
continuous basis from a heart tissue in contact with a catheter lead 510 to which the 
sensing circuit 520 is coupled, hi one embodiment, the sensed cardiac signal is an 

10 atrial signal In another embodiment, the sensed cardiac signal is a ventricular signal. 
Sensed cardiac signal from the sensing circuit 520 is then received and processed by 
an analyzer 550 of the controller 540 to compute curvatures contuiuously on a real- 
time basis at sample point Xi, X2, X3, Xi on the sensed cardiac signal. In some 

embodiments, the sample points are spaced at fixed intervals of time preset by the 

15 analyzer 550 on the sensed cardiac signal In one embodiment, the analyzer 550 

computes curvature at sample points Xi, X2, X3, Xi by fitting a cubic least 

square error curve to the sensed cardiac signal using N sample points to fit the cubic 
least square error curve. In one embodiment, the analyzer 550 computes the 
computing curvature (K) at a sample point Xi on the cardiac signal when using 5 

20 sample points over which the cubic fit is (N being 5 in this example) made using 

K = (2Ci/(l+B,y^) 

where 

Bi = TiY(I-2) + T2Y(M) + T3Y(I) + T4Y(I+1) + T5Y(I+2) 
25 Ci = SiY(I-2) + S2Y(I-2) + S3Y(I) + S4Y(H-1) + S5Y(I+2) 

where Si, S2, S3 ... and Ti, T2, T3 ...are constant terms in the summation. In this 
example embodiment, I is equal to the point at which the curvature is being computed. 
For example, if the curvature is being computed at sample point X3, then I will be 
equal to 3 and so on. In one embodiment, N is an odd number greater than or equal to 
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5. In this embodiment, the analyzer 550 computes curvature at a midpoint of the N 
number of sample points. 

In another embodiment, the analyzer 550 computes curvature based on an 
average curvature between two adjacent sample points based on a linear interpolation 
5 of the Bi and Ci between two adjacent points (integrating the computed curvatures 
between the two adjacent sample points). In one embodiment, the analyzer 550 
computes the curvatures continuously on a real-time basis from the sensed cardiac 
signal. Then the analyzer 550 extracts features by comparing the computed curvatures 
to a set of predetermined threshold values. In one embodiment, the predetermined 

10 threshold values are based on a previous curvature value, a curvature value, and a 
curvature threshold limit. In one embodiment, the analyzer 550 extracts features on a 
continuous basis. Then the analyzer 550 separates a set of features associated with a 
first complex upon detecting a second subsequent complex from the sensed cardiac 
signal. Then the analyzer 550 identifies a fiducial feature from the set of separated 

15 features associated with the first complex based on a predetermined deviation value. 
Then the analyzer 550 fixrther aligns the set of separated features with respect to the 
identified fiducial feature. In one embodiment, the analyzer 550 separates the set of 
features associated with the first complex based on a predetermined time earlier than 
the detected second subsequent complex. In one embodiment, the predetermined 

20 deviation value is based on a sample point having an amplitude farthest from a 
predetermined reference point. 

Then a comparator 570 coupled to the analyzer 550, compares the aligned set 
of features associated with the first complex with a predetermined set of templates, 
and classifies the first complex based on the outcome of the comparison. In one 

25 embodiment, the comparator 570 issues a command signal based on the outcome of 
the classification. In another embodiment, the therapy circuit 590 coupled to the 
comparator 570 delivers the electrical energy through the lead 510 to at least one 
electrode disposed in or around the heart 515 upon receiving the command signal from 
the comparator 570. The electrical energy delivered by the therapy circuit 590 can be a 

30 pacing pulse electrical energy. In one embodiment, the analyzer 550 includes a 
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variable gain circuitry that can adapt according to the changes in the sensed cardiac 
signal. In one embodiment, the controller 540 includes a memory 560 coupled to the 
comparator 570 and the analyzer 550 to store the extracted features of the first 
complex. In another embodiment, the memory 560 stores the classified first complexes 
5 for diagnostic purposes. 

It is to be understood that the above description is intended to be illustrative, 
and not restrictive. Many other embodiments will be apparent to those of skill in the 
art upon reviewing the above description. The scope of the invention should, therefore, 
be determined with reference to the appended claims, along with the full scope of 
10 equivalents to which such claims are entitled. 

Conclusion 

The above described method provides, among other things, a curvature based 
complex identification and classification. The process comprises sensing a cardiac 
15 signal and computing curvatures at sample points on the sensed cardiac signal, and 
then to extract features from the computed curvatures. Then to compare the extracted 
features with a set of predetermined templates, and then to classify the cardiac signal 
based on the outcome of the comparison. 

20 
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